Least-Squares Finite Element Methods for Quantum Electrodynamics

A significant amount of the computational time in large Monte Carlo simulations of lattice field theory is spent inverting the discrete Dirac operator. Unfortunately, traditional covariant finite difference discretizations of the Dirac operator present serious challenges for standard iterative methods. For interesting physical parameters, the discretized operator is large and ill-conditioned, and has random coefficients. More recently, adaptive algebraic multigrid (AMG) methods have been shown to be effective preconditioners for Wilson’s discretization [1] [2] of the Dirac equation. This paper presents an alternate discretization of the 2D Dirac operator of Quantum Electrodynamics (QED) based on least-squares finite elements. The discretization is systematically developed and physical properties of the resulting matrix system are discussed. Finally, numerical experiments are presented that demonstrate the effectiveness of adaptive smoothed aggregation (αSA ) multigrid as a preconditioner for the discrete field equations.